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Gas motions in the plane of the spiral galaxy NGC 3631 
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ABSTRACT 

The velocity field of the nearly face-on galaxy NGC 3631, derived from observations 
in the Ha line and Hi radio line, is analysed to study perturbations related to the 
spiral structure of the galaxy. We confirm our previous conclusion that the line-of- 
sight velocity field gives evidence of the wave nature of the observed two-armed spiral 
structure. Fourier analysis of the observed velocity field is used to determine the 
location of corotation of the spiral structure of this galaxy, and the radius of corotation 
R c is found to be about 42", or 3.2 kpc. The vector velocity field of the gas in the plane 
of the disc is restored, and taking into account that we previously investigated vertical 
motions, we now have a full 3D gaseous velocity field of the galaxy. We show clear 
evidence of the existence of two anticyclonic and four cyclonic vortices near corotation 
in a frame of reference rotating with the spiral pattern. The centres of the anticyclones 
lie between the observed spiral arms. The cyclones lie close to the observed spirals, 
but their centres are shifted from the maxima in brightness. 

Key words: galaxies: individual: NGC 3631 - galaxies: ISM - galaxies: kinematics 
and dynamics - galaxies: spiral - galaxies: structure - Ha line, 21cm line: galaxies. 
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1 INTRODUCTION 

The branch of astronomy known as dynamics of galactic 
discs has acquired, through long years of development, an 
unquestionably classical status, but has so far evolved as 
a part of theoretical astrophysics (e.g., Fridman & Poly- 
achenko 1984; Binney & Tremaine 1987). Observational 
data, which are the basis for dynamical investigations, have 
been up to now mostly one-dimensional: a mass distribu- 
tion in a disc is usually reconstructed from a surface bright- 
ness profile and from a long-slit major-axis velocity profile 
or rotation curve. Such an approach leads dynamicists to 
suppose a strict axisymmetry of galactic discs. However, ev- 
idently the latter are not axisymmetric in general. Particu- 
larly, bars and spiral arms are a clear manifestation of such 
non-axisymmetry. Thus, in order to make dynamical analy- 
ses more reliable, one needs two-dimensional (2D) data. 

With the advent of CCD detectors, 2D photometric 
studies began to appear. In particular, Kent (1984, 1985) 
has undertaken 2D decomposition of CCD images of galax- 



ies and has determined a lot of exponential disc parameters. 
In a series of works Athanassoula and co-workers (Considere 
& Athanassoula 1988; Garcia Gomez & Athanassoula 1993) 
have used an azimuthal Fourier analysis of images of galac- 
tic discs to reveal properties of their spiral structure, such as 
the number of arms and their pitch angle. Two-dimensional 
velocity fields, however, are rarely included into state-of- 
art dynamical investigations. The maximum yield obtained 
from such data is usually a rotation curve calculated in zero- 
order approximation of circular rotation, i.e. again under 
the assumption of axisymmetry. However, 2D velocity fields 
contain much more information. 

The observed line-of-sight velocity of gas in spiral galax- 
ies contains a contribution not only of the regular rotation, 
but also of the velocity perturbations due to the spiral den- 
sity waves. In principle, an analysis of the velocity field en- 
ables one to separate all these components of gas motion, but 
this task is far from simple: the expected amplitude of the 
main harmonics related to wave motion is about one order 
of magnitude lower than the maximal velocity of rotation of 
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a galaxy. In addition, the presence of both nonplanar oscilla- 
tions of the gas along the rotation axis and local non-circular 
motions makes the observed velocity field very complicated 
and difficult to interpret. Although spiral-related perturba- 
tions of the gas motion were detected both in our Galaxy 
(Yuan 1969 and references therein) and in many external 
galaxies beginning with the classical work of Rots on M81 
(1975), the amplitude of perturbed velocities and pattern 
angular velocity are badly known even for the best observed 
galaxies. 

Different methods were proposed to determine kinemat- 
ical parameters of density waves from the observed line-of- 
sight velocity fields (see Sakhibov & Smirnov 1987, 1989, 
1990; Bonnarel et al. 1988; Canzian, 1993; Sempere et al. 
1995; Schoenmakers et al. 1997; Westpfahl 1998 and refer- 
ences therein). However, all these methods have one or two 
principal shortcomings. First, they are based upon an al- 
leged possibility to restore the equilibrium rotation velocity 
without the analysis of the residual velocities. Such a possi- 
bility exists if we deal with the results of a model experiment 
and know the form of the gravitational potential (exactly 
the case with Canzian's 1993 investigation). However, when 
we analyse the line-of-sight velocity field of a real galaxy, 
independent restoration of the rotation curve becomes im- 
possible (Lyakhovich et al. 1997; Fridman et al. 1997). A 
rotation curve determined in the frame of a model of pure 
circular motion has systematic errors of the order of the 
residual velocities, and thus the residual field built on the 
basis of this curve does not represent the real field of veloc- 
ity perturbations caused by the density wave (Lyakhovich et 
al. 1997; Fridman et al. 1997). Second, all approaches men- 
tioned above are based on the assumption of a 2D character 
of the galactic motion in a disc, whereas any real galactic 
disc is a 3D object and regular motions induced by a den- 
sity wave are also three dimensional in principle (Fridman 
et al. 1997; Fridman et al., 1999). Thus the only direct ap- 
proach to analyse the observed velocity field is to seek for 
self-consistent solutions for the full vector velocity field. In 
other words, the rotation velocity and all three components 
of the perturbed velocity should be determined simultane- 
ously from the analysis of observational data, taking into 
account the 3D nature of the galactic discs. 

A recent attempt to restore the complete (three com- 
ponent) vector velocity field in the gaseous discs of grand- 
design galaxies from the observed field of line-of-sight ve- 
locities (Lyakhovich et al. 1997; Fridman et al. 1997), gave 
us a hope to build, in the future, an observationally-based 
foundation of the dynamics of the galactic discs. The knowl- 
edge of the complete velocity field gives us at once (1) rota- 
tion curve, (2) all the basic resonances: Lindblad and coro- 
tational, and (3) knowledge of the residual velocity field, 
containing recently discovered structures such as giant anti- 
cyclones (Fridman et al. 1997) and cyclones (Fridman et al. 
1999), the "constituent parts" of the spiral density waves. 
Finally, knowledge of the complete velocity field helps to 
determine the collective process — a kind of instability — 
which is responsible for the spiral-vortex structure of a given 
disc. Thus, we can state without risk of exaggeration that 
the observed velocity field provides the necessary observa- 
tional base for the construction of the dynamical portrait of 
a galaxy under consideration. 

The aim of the present article is to restore and anal- 



yse the velocity field of the gaseous disc of the grand-design 
galaxy NGC 3631, for which two types of line-of-sight veloc- 
ity data were obtained, well complementing each other: in 
the radio Hi and optical Ha lines (see below) . The Hi obser- 
vations used for this study were obtained by Knapen (1997) 
with the Westerbork Synthesis Radio Telescope, and the Ha 
observations were carried out at the Special Astrophysical 
Observatory (SAO) with its 6-m reflector equipped with an 
F/2.4 focal reducer and a scanning Fabry-Perot interferom- 
eter. 

NGC 3631 is a rather bright non-barred galaxy with 
well-defined spiral structure. Its optical axial ratio is close 
to unity: according to the RC3 catalogue (Vaucouleurs et al. 
1991), log a/b = 0.02 ± 0.07, so this galaxy looks nearly 
face-on. Such an orientation is very favourable for studying 
gas motions perpendicular to the plane of the galaxy, which 
was the main topic of our previous paper (Fridman et al., 
1998, referred hereafter as Paper I). 

In Paper I we showed that non-circular gas motions 
in NGC 3631 have a regular character, and that they are 
related to the observed two-armed spiral structure, which 
has a wave nature. The vertical (that is perpendicular to the 
plane of the disc) component of the gas motions as revealed 
by a Fourier analysis method (Fridman et al. 1997), was 
also found to be induced by the spiral density wave. The 
inclination angle of the disc of NGC 3631 was found to be 
about 17°, which enables, using the same observational data, 
the restoration of the vector velocity field in the plane of this 
galaxy, which is the main objective of the present paper. 

NGC 3631 is a grand-design spiral galaxy of type SAc, 
at a distance of 15.4 Mpc, as estimated from its recession 
velocity using a Hubble constant of 75 kms -1 Mpc -1 , which 
gives an angular scale of 75 pc per arcsec. Interestingly, the 
galaxy has been included in Arp's (1966) atlas of peculiar 
galaxies, thanks to its "straight arms" , and "absorption tube 
crossing from inside to outside of southern arm" . These fea- 
tures can be recognised in the i?-image shown in Fig. 1, but 
the galaxy as a whole looks to us rather normal. The atomic 
hydrogen distribution has been described by Knapen (1997 
and references therein to earlier work), and the ionized hy- 
drogen has been studied through emission in the Ha line by, 
among others, Boeshaar & Hodge (1977), Hodge (1982) and 
Rozas, Beckman & Knapen (1996). 

In Section 2 the results of a Fourier analysis of the 
observed distributions of optical (Ha and 7?-band) surface 
brightness and of Hi surface density are given, and compared 
with a Fourier analysis of azimuthal distributions of the ob- 
served line-of-sight velocities. Section 3 presents a model of 
large-scale gas motion, in which we assume that the gas ro- 
tates in the galactic plane, and simultaneously participates 
in the perturbed 3D-motions caused by the two-armed den- 
sity wave. The latter fact allows us to restrict the expansion 
of the line-of-sight velocity field to the first three Fourier 
harmonics (m obs = 1 — 3) (Fridman et al. 1997). Two in- 
dependent methods are used to determine the positions of 
corotation and other resonances, based on the relationships 
between the phases of azimuthal F v and radial F r oscil- 
lations of the perturbed velocity and the phase F a of the 
perturbed surface density. In Section 4 we then proceed to 
restore the vector velocity field of the gas in the plane of the 
disc. We show that in a frame of reference corotating with 
the spiral density wave pattern, giant cyclones exist along- 
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side anticyclones near the corotation radius. Anticyclones 
were described in our earlier work (Fridman et al. 1997), 
whereas cyclones were predicted (Fridman et al. 2000) in 
those gaseous galactic discs in which the gradient of the az- 
imuthal residual velocity exceeds the gradient of the rotation 
velocity in the reference frame corotating with spirals. We 
briefly summarize our main conclusions in Section 5. 

We refer to Paper I for a description of the observations 
which are used in this work. 



2 SPIRAL STRUCTURE OF NGC 3631 

Although the two-armed spiral structure of NGC 3631 is 
well defined in optical light, it has a number of irregulari- 
ties, especially in the outer parts of the disc. The Hi map 
has much lower angular resolution than optical images, yet 
it also clearly shows evidence of the two-armed spiral struc- 
ture (Knapen 1997). For the purpose of our study we should 
be convinced that the second harmonic of the brightness, 
mi, = 2, of the spiral structure exceeds all other harmon- 
ics. To check this, we divided the galactic disc into elliptical 
rings, corresponding to circular rings after deprojection, and 
carried out a Fourier analysis of the azimuthal brightness 
distribution, using an Ha image of NGC 3631 obtained at 
the SAO 6-m telescope through interferometric observations 
(Paper I), an i?-band image of the galaxy, as obtained from 
the Isaac Newton Group (ING) archive, and the neutral hy- 
drogen distribution, as obtained from 21 cm observations 
(Knapen 1997). We show the _R-band image, taken with the 
1-m Jacobus Kapteyn Telescope, in Fig. 1, which outlines 
the main spiral arm structure in this galaxy. Technical in- 
formation on this image can be found in Knapen (1997). 
The series of histograms in Fig. || shows the contributions of 
the different Fourier harmonics to the deviation from an axi- 
ally symmetrical distribution of brightness, or dispersion, for 
the three images mentioned above. Throughout the paper, 
we restrict our consideration to the region of the galactic 
disc R < 80" which corresponds to the extent of the optical 
spirals. 

Fig. ^ shows clearly that the second harmonic, which 
corresponds to the observed two-armed structure, indeed 
dominates the spectrum. The high level of the first harmonic 
in the Ha image is caused by the non-symmetrical distribu- 
tion of star-forming regions in the spiral arms, and does not 
reflect the true contribution of the first Fourier harmonic 
to the mass distribution in the galaxy. This assumption is 
supported by the low level of the first harmonic in the R- 
band image and the Hi map. To demonstrate this yet more 
clearly, we show in Fig. ^| the existence of a tight correlation 
between lines of maximum values of the second harmonic in 
the R- and Ha images (top panel), and a very weak corre- 
lation between points and lines of maximum values of the 
first harmonics (bottom panel). 

As shown earlier (Sakhibov & Smirnov 1989; Canzian 
1993; Fridman et al. 1997), if the circular velocity of gas in 
a galaxy is perturbed by a two-armed spiral pattern, this 
must lead to the appearance of the first and third Fourier 
harmonics (m b s = 1 and 3) in the azimuthal distribution 
of the observed line-of-sight velocity. In addition, the second 
harmonic (m b a = 2) may also appear if the density wave 
induces vertical oscillations of the gas (Fridman et al. 1997; 
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Figure 2. Contribution of individual Fourier harmonics to the 
deviation, or dispersion, of the brightness distribution from axial 
symmetry, as derived from (a) our H« image (Paper I), (b) our 
ii-band image (ING archive), and (c) 21 cm map (Knapen 1997). 



Fridman et al. 1998). The predominance of the first three 
harmonics in the line-of-sight velocity field of the galaxy is 
clear, as first demonstrated in paper I. Here, we show this 
result for both optical and radio velocity measurements in 
Fig. ^ and ^, using PA = 336° , which gives minimum dis- 
persion in a model of pure circular motion for the radio 
line-of-sight velocity data in the part of the galaxy under 
consideration. The optical data are practically insensitive to 
the change of PA by less than 10° . Fig. ^| shows the contri- 
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Figure 1. i?,-band image of NGC 3631 as obtained from the ING archive 



butions of different Fourier harmonics to the dispersion in 
the model of pure circular motion, averaged over the inner 
part of the disc (i? < 40"), whereas Fig. ^ shows similar 
histograms for the outer part of the optical disc (40" < R 
< 80"). 

Comparison of the Ha (Figs, ^a and |s]a) with the radio 
data (Figs, ^c and |H|c) shows that the amplitude of the 
harmonics is much smaller in the latter than in the former 
case. Such a difference could be a reason to distrust the 
optical or radio data, if these data were related to the same 
velocity field, but this is, in fact, not the case. 

As pointed out in Paper I, reasons for differences be- 
tween radio and optical estimates of the amplitudes of the 
second Fourier harmonics on the one hand, and of the first 
and the third harmonics on the other hand, should be differ- 
ent. In the latter case, the difference may be caused by the 
low resolution of the radio data. This is well illustrated in 
Figs. 4b and 5b, where the amplitudes were calculated after 
smoothing of the optical velocity field to a resolution of 14", 
close to that of the radio data. The squares of the ampli- 



tudes of the third harmonic in Fig. 4b and of the first and 
third harmonics in Fig. 5b are one-third/one- fifth as many 
as those in the original Ha data. It thus follows from Figs. 4 
and 5 that the histograms of the smoothed Ha data occupy 
an intermediate position between histograms of original Ha 
and Hi data — amplitudes in the histograms b) is closer to 
those in the histograms c) than those in histograms a). The 
second harmonic of the smoothed Ha line-of-sight velocity 
field is naturally much higher than the second harmonic of 
the Hi data. This can be explained by the different optical 
depth of the gaseous disc in the Ha and 21 cm lines, if one 
takes into account that the second harmonic is caused by 
the vertical motions in the density wave which are antisym- 
metrical with respect to the central plane of the disc (see 
discussion in Paper I). 

It follows from Figs. ^ and ^ that the contribution of 
harmonics with m Q b s > 3 into the observed line-of-sight ve- 
locity field is not significant for this galaxy. Note that this 
conclusion is not a universal rule - our preliminary analy- 
sis of data for other objects shows that it might be wrong 
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Figure 3. Superposition of the positions of the maximum value 
of the harmonics in the R (asterisks) and Ha (squares) images: 
(a) top panel, second harmonics; (b) lower panel, first harmonics. 



for some other galaxies, even though they have a two-armed 
grand-design spiral pattern. The anomalously large value of 
the amplitude of the 5th harmonic in Fig. is perhaps 
caused by the coincidence of the typical spatial scales of the 
5th harmonic in the outer part of the disc with the resolu- 
tion of the radio data. This hypothesis is supported by the 
appearance of an artificial anomalous amplitude of the fifth 
harmonic in the smoothed Ha data (Fig. pp). 

The extraction of the first three harmonics from the 
line-of-sight velocity field reveals that the residuals have a 
non-regular noisy-like character. Fig. ^| shows a histogram of 
the distribution of the velocity residuals with m > 3 all over 
the disc within the deprojected radius R < 80". It can be 
satisfactorily approximated by a Gaussian with a dispersion 
Vo^ of about 10 km/s, which is close to the mean pixel- 
to-pixel error in the individual velocity measurements. This 
allows the construction of a simplified model of the velocity 
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Figure 4. Contribution of individual Fourier harmonics of the 
line-of-sight velocity field into dispersion in the model of pure 
circular motion in the region R < 40", as derived from (a) original 
Ha line-of-sight velocity field, (b) Ha velocity field smoothed to a 
resolution comparable to the 21 cm data, (c) 21 cm velocity field. 



field of the galaxy, which uses only the first and the third 
harmonics to describe the gas motion in the plane of the 
disc. 

All results given below are restricted to this model. The 
second harmonic of the line-of-sight velocity field will be 
ignored here because it relates to the vertical gas motions 
and was considered earlier (see Paper I for details). 
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Figure 5. As Fig. 4, now for the region 40" < R < 



3 PHASE RELATIONSHIPS AND THE 
POSITION OF COROTATION 

The method used to restore the vector velocity field and 
to determine the corotation radius R c from observations of 
line-of-sight velocity distributions was described in detail 
by Lyakhovich et al. (1997) and applied to some galaxies 
by Fridman et al. (1997) and Fridman et al. (2000). The 
method is based on the comparison of Fourier coefficients 
of azimuthal distributions of the observed line-of-sight ve- 
locity, with those expected for a model where the perturbed 
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Figure 6. Histogram of residuals (differences between the ob- 
served and the model line-of-sight velocity fields) in the model of 
the line-of-sight velocity field which takes into account only the 
first three Fourier harmonics. The line shows the x _ distribution 
with a dispersion equal to that of the model, (a) For the origi- 
nal, unsmoothed, Ha line-of-sight velocity field, (b) for the 21 cm 
velocity field. 



velocity components at the particular moment of time V r , 
V v , V z are caused by a two-armed spiral density wave: 



V T (R, <p) = C r (R) cos[2tp - F r (R)] , 
V V (R, ip) = C 9 (R) cos[2cp - F^R)] , and 
V Z (R, ifi) = C Z (R) cos[2vj - F Z (R)] , 



(1) 
(2) 
(3) 



where d(R) and Fi(R) are an amplitude and phase of i- 
component of velocity. 

The main idea behind this is the following. The line-of- 
sight velocity is connected with the velocity components of 
the gas by the relationship (Lyakhovich et al. 1997; Fridman 
et al. 1997) 



V obs (R,ip) = V 3 + V V (R, <p) cost/? sin i + 
+ V r (R, ip) sin <p sini + V z (R, ip) cos i , 



(4) 



where V s is the systemic velocity of the galaxy. Taking into 
account that V r = V r , V v = V ro t + V v , V z = V z , where V ro t 
is the rotation velocity, and substituting (|l|) - (g) in (Q) we 
obtain the model representation of the line-of-sight velocity: 

mod ( 



V 



(R, <p) = V s + smi[ai(R) cosp + b\(R) sm tp + 
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+ a 2 (R) cos2v? + b 2 (R) sin2<p + 
+ 03 (R) cos3v? + 63 (R) sin3(y3] , 



(5) 



where Fourier coefficients related to phases and amplitudes 
of the velocity components are: 

C r sin F r + Gin cos Fin 



ffll = Vrot + 



61 = - 



C r COS Fr 



2 

C v sin F w 



«2 
&2 

O3 



C z cos F z cot i , 

C z sin F z cot i , 

C r sin F r — Cp cos F v 
2 

C r cos FV + Cm sin F w 



(6) 

(7) 

(8) 
(9) 

(10) 
(11) 



Calculating the Fourier coefficients of the observed line- 
of-sight velocity field (a° ba , b° ba , with i = 1, 2, 3) and equat- 
ing them to the model ones (Eqs. p[|i"l|), we obtain the base 
to determine the amplitudes and phases of all three velocity 
components. 



The Fourier coefficients a° bs , b° b ° , af" , b? 3 , a?" and 



63 s as well as the best fit parameters of the galactic disc 
(PA, inclination, and centre position) may be found from 
the observed line-of-sight velocity distribution V° bs (r j , ip j) 
by minimising in each elliptical ring the quantity X 2 (Rk), as 
determined by the following equation (see Lyakhovich et al. 
1997; Fridman et al. 1997 for details): 

x 2 (R k ) = J2 <yi( R i><Pi) - v ° - 
3 

n=3 \ 2 

— S ' [a° bs (R k ) cos nip j + b^ s (R k ) sinnp^ sini ] ,(12) 

n=l / 

where k designates the number of the ring and summing is 
performed over all pixels belonging to the ring. 

The first harmonic of the observed velocity field con- 
tains contribution from both the rotation velocity and the 
perturbed motion, which can not be separated without tak- 
ing some additional proposition, but the third harmonic 
should be unambiguously related to the velocity perturba- 
tion connected with the observed spiral arms. To verify this, 
Fridman et al. (1997) proposed to use the "modified third 
harmonic", which has a form cos(2y> — F3 + 7r/2), where F3 
is the phase of the original third harmonic of the observed 
line-of-sight velocity. In the case of a tightly wound spiral, 
it was shown that the maxima of this "modified third har- 
monic" follow the maxima of the perturbed surface density 
of the disc outside corotation (that is F3 = F a -+- 7r/2). In- 
side corotation, they can coincide with maxima or minima 
of the perturbed surface density, depending on the relation 
between the amplitudes of the radial and azimuthal residual 
velocities. If CV < C v then F3 = F a + 7r/2, otherwise F3 = 
F(j — 7r/2. Fig. [?] shows the correlation between the positions 
of maxima of the perturbed surface density (characterized 
by the m b s = 2 harmonic in the i?-band brightness dis- 
tribution) and maxima (a) and minima (b) of the "modified 
third harmonic". As one can see, the latter set of maxima 
really follow the spiral arms in the outer part of the disc 
and in the very inner region, whereas they lie between the 
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Figure 7. Comparison of the positions of the maxima (a) and 
minima (b) of the "modified third harmonic" of the Ha line-of- 
sight velocity field (triangles) with the form of the spiral arms 
characterized by positions of maxima of the second Fourier har- 
monic of R brightness map (asterisk). The modified third har- 
monic is cos(2i/3 — F3 + 7r/2), where F3 is the phase of the original 
third harmonic. 



arms in the radial region 25" < R < 40". This proves the 
correctness of the assumption we used on the wave nature 
of the spiral structure in NGC 3631, and implies that the 
position of corotation is at about 40". This argumentation 
gives the first rough estimate of the corotation radius R c . 

To obtain a more accurate estimate of the position of 
corotation we previously proposed two methods (Lyakhovich 
et al. 1997; Fridman et al., 1997). The first method is based 
on the comparison of the radial behaviour of the sinus com- 
ponents of the first (b° bs ) and third (b^ bs ) Fourier harmonics 
of the line-of-sight velocity field. In the case of tightly wound 
spirals, the following relations are fulfilled (for details see 
Lyakhovich et al. 1997 and Fridman et al., 1997): 



© RAS, MNRAS 000, 



8 A.M.Fridman et al. 
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c) 




Figure 8. Behaviour of \b° b3 (R)\ - \b° bs (R)\ as a function of 
galactocentric radius R in NGC 3631, as derived from (a) origi- 
nal Ha line-of-sight velocity field, (b) smoothed Ha velocity field, 
(c) 21 cm velocity field. Error bars correspond to 3cr level. Ac- 
cording to results of the density wave theory in the approxima- 
tion of tightly wound spirals, the difference should be negative 
inside corotation and positive outside. Thus the data show that 
the corotation radius is at 40" ± 7". 



'(R)\-W bs (R)\<0, 



\b? s (R)\ 



\b? 3 (R)\>0, 



for 
for 



R < R c , and 
R> R c - 



(13) 



These inequalities enable the determination of the lo- 
cation of the corotation radius from observational data. Ac- 
cording to Eq.([L3|), corotation is located in the region where 
the difference |6f "{R)\ - \b° b S {R)\ chan ges sign, from minus 
to plus. 
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Figure 9. Variation of (bSj 68 
dius R in NGC 3631, as derived from (a) original Ha line-of-sight 
velocity field, (b) smoothed Ha velocity field, (c) 21 cm velocity 
field. Error bars correspond to 3cr level. According to results of 
the density wave theory in the approximation of tightly wound 
spirals, the difference is negative inside and positive outside coro- 
tation. Thus the data show the corotation radius to be about 43" 
± 4". 



- b° bs ) cos F a with galactocentric ra- 



Fig. |shows the radial behaviour of |&3 bs (ii)| - \b° b S (R)\ 
in NGC 3631, as obtained from the original Ha (a), 
smoothed Ha (b), and 21-cm (c) line-of-sight velocity fields. 
Fig. |^ shows that this function derived from Ha line-of-sight 
velocity field is negative within the errors in the inner part 
of the galaxy, and positive in the outer region, in accordance 
with the expectations from ([H5]). From these data it follows 
that the corotation radius is at R — 40" ± 7". The level of 
errors in 21-cm data is relatively higher, due to lower resolu- 
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tion. Nevertheless, behaviour of \b% S {R)\ 



|&f s (7?)| derived 



from 21-cm line-of-sight velocity is in good agreement with 
Ha data within one a level of confidence. 

In addition, yet another method to estimate R c was pro- 
posed (Lyakhovich et al. 1997; Fridman et al. 1997), based 
on the relation between the phases of the perturbed surface 
density and the radial perturbed velocity of a gas, which is 
fulfilled for trailing tightly wound spirals: 



Fa F r = 



7r at R < R c 
at R > R c 



= (1 — signtD) 7r/2 . 



(14) 



This relation shows that gas should move along spiral 
arms inwards inside of corotation, and outwards outside of 
it. As shown in the cited papers, it may also be written as 
the relationship between the phase of the perturbed surface 

and6f s : 



density and the Fourier coefficients 63 
(bf s (R) -bf 3 (R)) cos F a {R) < 0, 



(R)-bT s (R)) cos F„(R) 



for R< R c 
for R > R c 



(15) 



Fig. I shows the dependence of {bf s - bf 3 ) cos F a on 
galactocentric radius R. The situation is quite similar to 
presented in the previous figure. The current approach leads 
to the estimate of R c of 43" ± 4". 

Combining the results from these methods, which agree 
well, we may conclude that the corotation radius R c ~ 42" ± 
5", or 3.2 kpc±0.38 kpc. 



4 RESTORED VELOCITY FIELD OF GAS IN 
THE GALACTIC PLANE 

As shown in Section 2, lowering the angular resolution leads 
to an underestimate of the amplitudes of the Fourier har- 
monics. Therefore, to restore a two component vector veloc- 
ity field of a gas in the disc of NGC 3631, we use the full 
resolution Ha line-of-sight velocity field. 

To restore the velocity field in the plane of a galactic 
disc, i.e., its radial V r and azimuthal V v = Vrot + Vy, velocity 
components, it is necessary to determine five unknown func- 
tions: Vrot(r), C r (r), C v (r), F r (r), F v (r) (see Eqs. | and|). 
These five functions are connected with the Fourier coeffi- 
cients of the observed line-of-sight velocity field by the four 
relations (^), ([?]), (|lo|), and (p"i"|). An additional condition, re- 
quired to close the system, should have a theoretical origin. 
Unfortunately, up to now, a reliable condition, valid for any 
density wave amplitude, is not available. Several possibilities 
discussed in the literature (Sakhibov & Smirnov, 1989; Frid- 
man et al. 1997) are based on some kind of approximation 
and have limited applicability. To overcome this difficulty, 
we propose the following approach. 

Among the functions listed above, V ro t(r) can be most 
reliably estimated in two independent ways from observa- 
tional data. The first one uses the equilibrium condition of 
a gaseous disc rotating in a gravitational potential 

V r 2 ot /r = d*/dr. (16) 

The right-hand side of the Eq.(jl|) is determined from 
the mass distribution in a galaxy, or its surface bright- 
ness maps, assuming the mass-to-light ratio is known and 
constant with radius. For this purpose, we use a three- 
component dynamical model of a spiral galaxy similar to the 
one used by Sumin, Fridman, & Haud (1991). Although the 
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Figure 10. The radial dependence of b^ bs — b° bB = C' r cos 
b obs + b obs = c ^ sin observed in the spiral galaxy NGC 3631. 
An estimate of the amplitudes of the velocity components from 
the extremes of these functions gives max(C r ) — max(C v ) ~ 
60 km/s. 



model is rather crude, the resulting rotation curve V ro t(R) 
corresponds to the region approximately between (ai) m i n 
and (ai) max . 

The same result can be obtained in another way. From 
equation (^|) it follows that the difference \ai — V ro t\ can- 
not exceed the amplitudes CV and C v .which, in turn, are 
connected by equations (§), (£(]), and (y) with the Fourier 
coefficients bi bs ^a3 bs , and b^ bs , determined from our obser- 
vations. In Fig. hoJ, we show the radial behaviour of b^ bs — b° bs 
= C r cos F r and bf s + b° ba = sin F v . The extremes of 
these functions allow estimates of the amplitudes C- and 
CV From Fig. we conclude that in NGC 3631 a maxi- 
mum value of the amplitude of the residual velocities occurs 
at 60 km/s, i.e. 

|ai - V ro t\max < 60 km/s. (17) 

The conditions ( |l6| ) and ([n]) do not allow an exact cal- 
culation of the function V ro t{r). Nevertheless, they set lim- 
its on the variations of both the amplitude and the form of 
V ro t(r). Within these limits, varying the mass-to- light ratio 
(±40%) and assuming it does not depending on radius we 
obtain a set of trial curves (Fig. |ll|) and analyse the velocity 
fields restored from Eqs. (|E|)-(|ll|) for a given V ro t(r). 

Fig. [l^ shows four examples of the restored velocity 
field of NGC 3631, in the plane of the disc, and in a ref- 
erence frame rotating with the spiral pattern. Positions of 
the maxima of the second Fourier harmonic of the _R-band 
surface brightness map are overlaid to indicate the relative 
location of spiral arms and structures in the velocity field. 

It is clearly seen that, in spite of broad variations in 
Vrot{r), the general structure of the velocity field changes 
only slightly. In all cases this structure demonstrates the 
presence of two anticyclones, located on the corotation cir- 
cle, and between spiral arms. Variations of the rotation curve 
only lead to small changes in the basic quantitative features 
of the anticyclones. In all cases, the maximum noncircular 
velocity in the vortices is about 60 km/s, their radial width 
is about 20", and the variation of the azimuthal position of 
the centres of anticyclones is less than 10° . This proves that 
the exact shape of the rotation curve does not impact on 
the basic conclusion of the existence of giant anticyclones in 
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Figure 12. Restored velocity field of NGC 3631 in the plane of the disc in the reference frame rotating with the pattern speed. Overlaid 
asterisks show the locations of maxima of the second Fourier harmonic of the ij-band brightness map of the galaxy. The circle marks 
the position of the corotation. Solid lines demonstrate vortex separatrices or nearly closed streamlines in the absence of a separatrix (the 
ones greater in size correspond to anticyclones, and the smaller - to cyclones), 
(a) The curve marked by symbol "a" in Fig. 11 used as V ro t(r). 



the gaseous disc of NGC 3631. These anticyclones are simi- 
lar to those revealed earlier in the velocity field of NGC 157 
(Fridman et al. 1997). 

Another new result is that regions of cyclonic shear are 
seen in Fig. |l2| As a consequence of the relatively high am- 
plitude of the density wave in NGC 3631, a cyclonic shear is 
produced by the density wave in some regions, which domi- 
nates the anticyclonic shear caused by differential rotation. 
The appearance of cyclones in gaseous galactic discs with a 
strong density wave was predicted earlier by Fridman et al. 
(2000). 



To choose between the velocity fields presented in 
Fig. ^ the field which is closest to the real velocity field 
of the galactic disk of NGC 3631, we use the following cri- 
terion. In the course of many revolutions the growth of the 
density wave amplitude is stopped by the saturation of a 
corresponding instability. Thus in the reference frame coro- 
tating with spiral arms the velocity field of the galactic disk 
should be stationary. In such a field the vortices should have 
closed separatrices (lines dividing two families of trajecto- 
ries: trapped and transit ones). According to the criterion 
stated above, the velocity field presented in Fig. [I2|c is dis- 
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Figure 12. (b) The curve marked by symbol "b" in Fig. 11 is used as V ro t{i')- 



tinct in that sense that it is only case where cyclones have 
close separatrices. It is an argument in favour of choosing 
just this example as the closest to the real velocity field of 
the NGC 3631 gaseous galactic disk. At the same time that 
means that the curve "c" in Fig. [H] is close to the real rota- 
tion curve of the disk. The relative position of the vortices 
and spiral arms in Fig. [L2]c is in a good agreement with 
theoretical predictions (Fridman et al. 1999). 



5 CONCLUSIONS 

We can briefly summarize our main conclusions from this 
paper as follows: 



(i) Based on an analysis of line-of-sight velocity fields of 
gaseous emission lines for the spiral galaxy NGC 3631 we 
confirm observationally the theoretical conclusions from our 
previous work about the wave nature of its two-armed spiral 
structure, and discuss its origins and properties in light of 
this theoretical framework. 

(ii) Using two independent methods, we find that the 
corotation radius in this galaxy is at about 42" or 3.2 kpc. 

(iii) The projection of the restored three-dimensional vec- 
tor velocity field of gas in the plane of the galaxy, and in a 
reference frame corotating with the spiral pattern, reveals 
the presence of two anticyclonic vortices near corotation. 
We thus confirm the theoretical prediction foreseen earlier 
on the basis of a study of the general principles of the wave 
nature of spiral structure in galaxies. 
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Figure 12. (c) The curve marked by symbol "c" in Fig. 11 is used as V ro t(r)- 



(iv) We show the existence of cyclonic vortices in 
NGC 3631, apart from that of anticyclones, as mentioned 
above. Such cyclonic vortices are a consequence of a high 
amplitude of the density wave in this galaxy. 
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Figure 11. Examples of trial curves used to represent the ro- 
tation curve (V ro t(r)) in NGC 3631, shown by solid lines, along 
with the observed behaviour of a° bs (r) (triangles). Abscissa is the 
radius in arcsec and ordinate — rotation velocity in km/s. Ac- 
cording to the analysis presented below, the thickest line marks 
the curve corresponding to the real rotation curve of NGC 3631 
gaseous disk. 
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